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5 INCORPORATION BY REFERENCE OF MATERIAL SUBMITTED ON 

COMPACT DISC 

A computer program listing appendix that includes a genetic algorithm utilized in 
accordance with aspects of the presently disclosed technology is contained on a submitted 
compact disc. Each of two identical copies of such compact disc includes a file named 
10 "CXU-33.9 Genetic Algorithm", dated 04/04/2005 and having a size of 126 KB. The 
program listing contained in such file is hereby incorporated by reference for all 
purposes. 

FIELD OF THE INVENTION 

15 This technology provides a method (application) of an algorithm to facilitate the 

design of wideband operations of antennas, and the design of sleeve cage monopole and 
sleeve helix, units. The technology is of interest/commercial potential throughout the 
audio communications community. 

Omnidirectional capabilities and enhanced wideband capabilities are two 

20 desirable features for the design of many antenna applications. Designing 

omnidirectional antennas with wideband capabilities requires rapid resolution of complex 
relationship among antenna components to yield an optimal system. The invention 
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comprises the use of a genetic algorithm with fitness values for design factors expressed 
in terms to yield optimum combinations of at least two types of antennas. 

Cage antennas are optimized via a genetic algorithm (GA) for operation over a 
wide band with low voltage standing wave ratio (VSWR). Numerical results are 
5 compared to those of other dual band and broadband antennas from the literature. 
Measured results for one cage antenna are presented. 

Genetic algorithms and an integral equation solver are employed to determine the 
position and lengths of parasitic wires around a cage antenna in order to minimize voltage 
standing wave ratio (VSWR) over a band. The cage is replaced by a normal mode 
10 quadrifilar helix for height reduction and the parasites are re-optimized. Measurements 
of the input characteristics of these optimized structures are presented along with data 
obtained from solving the electric field integral equation. 

Genetic algorithms (Y. Rahmat-Samii and E. Michielssen, Electromagnetic 
Optimizations by Genetic Algorithms, New York: John Wiley and Sons, Inc., 1999) are 
15 used here in conjunction with an integral equation solution technique to determine the 
placement of the parasitic wires around a driven cage. The cage may be replaced by a 
quadrifilar helix operating in the normal mode in order to shorten the antenna. 
Measurements of these optimized structures are included for verification of the 
bandwidth improvements. 

20 

BACKGROUND OF THE INVENTION 

Recent advances in modern mobile communication systems, especially those 
whichemploy spread-spectrum techniques such as frequency hopping, require antennas 
whichhave omnidirectional radiation characteristics, are of low profile, and can be 

25 operatedover a very wide frequency range. The simple whip and the helical antenna 
operating inits normal mode appear to be attractive for this application because they 
naturally haveomnidirectional characteristics and are mechanically simple. However, 
these structuresare inherently narrow band and fall short of needs in this regard. Hence, 
additionalinvestigations must be undertaken to develop methods to meet the wide 

30 bandwidthrequirement of the communication systems. 
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BRIEF SUMMARY OF THE INVENTION 

This invention comprises a method to design (produce) a product and the 
product(s) designed/produced as a result of the application of the method. The products 
are broadband, omnidimensional communications antennas, and the design procedure 
5 involves the coordinated, sequential application of two algorithms: a generally described 
"genetic algorithm that simulates population response to selection and a new algorithm 
that is a fast wire integral equation solver that generates optimal multiple antenna designs 
from ranges of data that limit the end product. Individual designs comprise a population 
of designs upon which specified selection by limiting the genetic algorithm ultimate 

10 identifies the optimum design(s) for specified conditions. Superior designs so identified 
can be regrouped and a new population of designs generated for further 
selection/refinement. 

The products are the antenna designs and specifications derived as a product of 
the application of the method briefly described above. The antennas all are characterized 

15 generally as broad band and omni directional, two features of critical importance in 
antenna design. In addition, although much of the theory has been developed on 
monopole antennas, both the method and designs include both monopole and dipole 
designs. In addition, the designs include sleeve-cage and sleeve-helix designs as 
hereinbelow further described. 

20 The cage monopole comprises four vertical, straight wires connected in parallel 

and driven from a common stalk at the ground plane. The parallel straight wires are 
joined by crosses made of brass (or other conductive) strips, the width of which is equal 
to the electrical equivalent of the wire radius. Compared to a single wire, this cage 
structure has a lower peak voltage standing wave ratio (VS WR) over the band. A 

25 structure with lower VSWR is amenable to improved bandwidth characteristics with the 
addition of parasitic elements. 

Adding parasitic elements of equal height and distance from the center of the cage 
monopole creates a sleeve cage monopole. The sleeve cage monopole has a greater 
bandwidth than its otherwise comparable antennas. Fitness values are determined by 

30 relative bandwidth, with greater fitness being associated with wider bandwidth defined by 
£2/fl, where f2 and fl are respectively the largest and smallest frequencies between 
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which VSWR is 3.5 or less. Speed of optimization is increased by interpolation of the 
impedance matrix. 

The heart of the process is the solution of the equation governing total axial 
current. The executable algorithm linked to the genetic algorithm by the fast wire integral 
5 equation solver provides a rapid method of solving this equation for varied values and 
inputs. The basic theory and equations are incorporated completely herein. See, S.D. 
Rogers and CM. Butler, "An efficient curved-wire integral equation solution technique," 
submitted to IEEE Trans. Antennas Propagat 

Reduced height without loss of bandwidth or omni-directional capabilities is a 

10 desired feature of antenna designs for a plurality of applications. These include 

installations in vehicles and confined interior spaces. The helix structure yields shorter 
antennas than the traditional whip structure with otherwise comparable features. Height is 
a function of the pitch angle of the helix, such that a pitch angle of 42 degrees reduces 
height by 30 percent. The addition of parasitic elements reduces VSWR in a helix 

15 configuration in a magnitude similar to the reduction noted for the cage monopole design. 

Many modern wireless communication systems require low-profile antennas. To 
meet this requirement, we consider the helical antenna operating in the normal mode. A 
normal mode helix and a straight wire antenna having approximately the same wire 
length exhibit similar input impedance and far field patterns. One drawback to the helix 

20 operating in the normal mode is that its bandwidth is too limited for many applications. 
To increase the bandwidth, we have considered several potential remedies, one of which 
is discussed herein. It is well known that adding additional parasitic straight wires 
on either side of a driven dipole antenna may increase the bandwidth of the dipole (J.L. 
Wong and H.E. King, AP-21, no. 5, 725-727, Sept. 1973). One must be especially 

25 careful, however, to choose parasitic elements with the proper length and spacing. We 
use this basic idea to increase the bandwidth of the helical monopole. A structure similar 
to the sleeve dipole but applied to the helix has been used to design dual frequency 
antennas (P. Eratuuli, et. al., Electronics Letters, v. 32, no. 12, 1051-1052, June 1996). 
The helix and its helical sleeve are both driven in the antenna of this reference. Several 

30 novel antenna structures are considered such as a driven helical antenna adjacent to 

parasitic helices and straight wires. Another candidate structure consists of a driven helix 
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with helical parasites inside or outside of the driven element, which has the added benefit 
of conserving space. In any case, due to the large number of parameters in a helix, it is 
more difficult to design a broadband sleeve helical antenna than is the case for a sleeve 
dipole. It is not feasible to obtain optimum values of parameters by trial and error. Thus 
5 we employ a genetic algorithm routine (D.L. Carroll, A FORTRAN Genetic Algorithm 
Driver, http://www.staff.uiuc.edu/-carroll/ga.html) and efficient integral equation 
solution techniques to optimize the antenna system for bandwidth. Having an efficient 
numerical solution technique is necessary for this problem since the geometry of the 
antenna is redefined for each structure evaluated by the genetic algorithm. Since these 

10 antennas have a high degree of curvature, their solutions generally require a large number 
of unknowns for representing the geometry. An efficient solution technique which gets 
around this problem is used (S.D. Rogers and CM. Butler, APS Symposium Digest, vol. 
I, 68-71, July 1997). 
1) Commission B. B-2 Antennas 

15 2) A genetic algorithm is used to optimize helical parasitic elements for a helical 
antenna. 

3) This work is an extension to increasing the bandwidth of dipoles by use of parasites. 
This research could not have been completed in a time efficient manner without the 
development of an efficient integral equation solution technique for curved wires with 
20 reference below. 



S.D. Rogers and CM. Butler, "Reduced Rank Matrices for Curved Wire 
Structures," Digest of IEEE APS Symposium, Montreal, Canada, vol. I, pp. 68-71, July 
1997. 

25 

We have recently shown from numerical calculations that the bandwidth of a 
normal mode helix can be increased by the addition of close-by wire parasites (S.D. 
Rogers, J.C Young, and CM. Butler, "Bandwidth Enhanced Normal Mode Helical 
Antennas," Digest 1998 USNC / URSI National Radio Science Meeting, Atlanta, GA., p. 
30 293, June 1998). Agenetic algorithm and a fast integral equation solution technique are 
employed to determine the optimum distance and height of these parasites. In (H.E. King 
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and J.L. Wong, "An Experimental Study of a Balun-Fed Open-Sleeve Dipole in Front of 
a Metallic Reflector," IEEE Trans. Antennas Propagat. (Commun.), vol AP-20, pp. 201- 
204, March 1972) these parameters were determined experimentally when the driven 
element was a straight wire dipole. In a recent paper (H. Nakano, et. al., "Realization of 
5 Dual-Frequency and Wide-Band VSWR Performances Using Normal-Mode Helical and 
Inverted-F Antennas," IEEE Trans. Antennas Propagat., vol. AP-46, June 1998) a central 
parasitic straight cylinder was added inside a driven single wire helix to obtain dual 
frequency operation. We have found that even greater bandwidth, over that of a single 
driven wire, can be realized when the parasites are placed around "cage" monopoles 

10 having several parallel wires. Similar observations are made about a multifilament 
versus a single filament helix. 

Bandwidth of vertically polarized wire antennas is often increased by adjustment 
of the antenna geometry. King and Wong reduce VSWR by placing parasitic wires 
around a driven element, creating the well-known open sleeve dipole (KING, H.E., and 

15 WONG, J.L.: 'An experimental study of a balun-fed open-sleeve dipole in front of a 
metallic reflector 5 , IEEE Trans. Antennas Propagat, March 1972, 20, (2), pp. 201-204). 
In (NAKANO, H., IKED A, N., WU, Y., SUZUKI, R., MIMAKI, H., and YAMAUCHI, 
J.: 'Realization of dual-frequency and wide-band VSWR performances using normal- 
mode helical and inverted-F antennas', IEEE Trans. Antennas Propagat., June 1998, 46, 

20 (6), pp. 788-793) the displacement of a parasitic monopole is varied inside a driven 

normal mode helical antenna in order to control its characteristics. Cage antennas can be 
made broadband when their dimensions are chosen judiciously. Genetic algorithms and 
integral equation solution techniques are employed here to optimize the dimensions of the 
cage antenna in order to create a structure with low VSWR over a wide band. 

25 Recent advances in modern mobile communication systems, especially those 

which employ spread-spectrum techniques such as frequency hopping, require low- 
profile, broadband, omnidirectional (in azimuth) antennas. The simple whip and the 
helical antenna, operating in its normal mode, are potentially attractive for these 
applications because they naturally have suitable radiation characteristics and are 

30 mechanically simple and rugged. However, these structures are inherently narrow band. 
Hence, additional measures are employed to meet the wide bandwidth requirement of 
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communication systems. Antennas often are loaded with tuning circuits and are 
connected to radios through matching networks in order to improve overall bandwidth. 
Altering the antenna geometry is another method for modifying bandwidth properties. 
The sleeve monopole, in which the outer conductor of the coaxial feed line forms a 
5 "sleeve" around the base of the protruding center conductor, is known to have greater 
bandwidth than the conventional monopole and has been studied extensively (J. Taylor, 
"The sleeve antenna," doctoral dissertation, Cruft Lab., Harvard Univ., Cambridge, MA, 
1950); (R.W.P. King, The Theory of Linear Antennas. Cambridge, MA: Harvard Univ 
Press, 1956); (A. J. Poggio and P.E. Mayes, "Pattern bandwidth optimization of the sleeve 

10 monopole antenna," IEEE Trans. Antennas Propagat. (Commun.), vol. AP-14, pp. 643- 
645, Sept. 1966); (Z. Shen and R. MacPhie, "Rigorous evaluation of the input impedance 
of a sleeve monopole by modal-expansion method," IEEE Trans. Antennas Propagate 
vol. AP-44, pp. 1584-1591, Dec. 1996). A variation of this antenna is the open-sleeve 
dipole which has straight- wire parasites in place of the coaxial sleeve. The effects of the 

15 spacing and size of the parasitic elements on the VSWR are determined experimentally in 
(H.E. King and J.L. Wong, "An experimental study of a balun-fed open-sleeve dipole in 
front of a metallic reflector," IEEE Trans. Antennas Propagat. (Commun.), vol. AP-20, 
pp. 201-204, March 1972). In other papers, parasitic and driven elements of various sorts 
are combined in order to create dual band antennas. In (P. Eratuuli, et. al. 9 "Dual 

20 frequency wire antennas," Electronics Letters, vol. 32, no. 12, pp. 1051-1052, June 6, 

1996) the driven wire is a straight monopole or a helix surrounded by a parasitic helix. In 
(H. Nakano, et. ah, "Realization of dual-frequency and wide-band VSWR performances 
using normal-mode helical and inverted-F antennas," IEEE Trans. Antennas Propagate 
vol. AP-46, pp. 788-793, June 1998) the position of a straight-wire parasite inside a 

25 driven normal mode helical antenna is adjusted to control the VSWR over the band of 
operation. Another antenna, which can be made to have broadband properties if its 
dimensions are chosen judiciously, is the cage antenna (S.D. Rogers and CM. Butler, 
"Cage antennas optimized for bandwidth," submitted to Electronics Letters, April 2000). 
The cage is more amenable than a single straight wire to improvement in bandwidth 

30 when parasitic wires of appropriate size and spacing are added (S.D. Rogers and CM. 
Butler, "The sleeve-cage monopole and sleeve helix for wideband operation," Digest of 
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APS Symposium, Orlando Florida, vol. 2, pp. 1308-1311, July 1999). 

We have found that the cage structure and multifilar helices are more amenable 
than single wire antennas to improvements in VSWR when parasitic wires are added. The 
helical configuration can be used to reduce the height of the antenna, but at the sacrifice 
5 of bandwidth. While the addition of the parasitic wires improves the overall bandwidth, 
the VSWR increases outside the design band. Fast integral equation solution techniques 
and optimization methods have been developed in the course of this work and have led to 
effective tools for designing broadband antennas. 

Certain exemplary attributes of the invention may relate to a method to create 
10 optimum design specifications for omni-directional, wide band antennas comprising the 
steps of: 

(a) loading software including a genetic algorithm and an executable algorithm that is a 
fast wire equation solver into a computer; 

(b) loading instructions into said computer specifying basic antenna design to be 
15 optimized; 

(c) loading antenna design parameters and corresponding ranges of values for said 
parameters into said computer; 

(d) specifying resolution of said parameters by loading number of bits per parameter into 
said computer; 

20 (e) executing (operating) said genetic algorithm thereby generating a population of 
individual antenna designs each with a fitness value; and 

(f) evaluating relative fitness of antenna designs produced and selecting superior designs 
for continued refinement. 

The foregoing method may further comprise the following exemplary subroutines 
25 and algorithms for the software involved: 

(a) a first algorithm that allows different values for critical design elements to combined 
in all possible combinations and a fitness value for each design ultimately estimated; 

(b) a second algorithm that determines electronic current in an antenna by solving an 
integral equation numerically; 

30 (c) a computer program link that provides essential communication between said first 
algorithm and said second algorithm. 
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Certain exemplary attributes of the invention may further relate to the sleeve 
monopole antenna designs, the cage sleeve monopole antenna designs, and the sleeve 
dipole antenna designs produced following the foregoing methods. Those of ordinary 
skill in the art will appreciate that various modifications and variations may be practiced 
5 in particular embodiments of the subject invention in keeping with the broader principles 
of the invention disclosed herein. The disclosures of all the citations herein referenced are 
fully incorporated by reference to this disclosure. 

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWINGS 

10 A full and enabling description of the presently disclosed subject matter, 

including the best mode thereof, directed to one of ordinary skill in the art, is set forth in 
the specification, which makes reference to the appended figures, in which: 

FIG. 1 A is a perspective view of an exemplary cage monopole antenna 
embodiment in accordance with the presently disclosed technology; 
15 FIG. IB is a graphical representation of the voltage standing wave ratio (VSWR) 

versus frequency for the exemplary cage monopole antenna of FIG. 1 A; 

FIG. 1C is a graphical representation of the directivity versus frequency for the 
exemplary cage monopole antenna of FIG. 1A; 

FIG. 2 A is a perspective view of an exemplary sleeve-cage monopole antenna 
20 embodiment in accordance with the presently disclosed technology; 

FIG. 2B is a graphical representation of the voltage standing wave ratio (VSWR) 
versus frequency for the exemplary sleeve-cage monopole antenna of FIG. 2A; 

FIG. 2C is a graphical representation of the directivity versus frequency for the 
exemplary sleeve-cage monopole antenna of FIG. 2 A; 
25 FIG. 3 A is a perspective view of an exemplary quadrifilar helix antenna 

embodiment in accordance with the presently disclosed technology; 

FIG. 3B is a graphical representation of the voltage standing wave ratio (VSWR) 
versus frequency for the exemplary quadrifilar helix antenna of FIG. 3 A; 

FIG. 3C is a graphical representation of the directivity versus frequency for the 
30 exemplary quadrifilar helix antenna of FIG. 3 A; 



FIG. 4A is a perspective view of an exemplary sleeve helix antenna embodiment 
in accordance with the presently disclosed technology;FIG. 4B is a graphical 
representation of the voltage standing wave ratio (VSWR) versus frequency for the 
exemplary sleeve helix antenna of FIG. 4A; 
5 FIG. 4C is a graphical representation of the directivity versus frequency for the 

exemplary sleeve helix antenna of FIG. 4 A; 

FIG. 5 A is a graphical representation of VSWR versus frequency for a cage 
antenna optimized for VSWR < 2.0; 

FIG. 5B is a graphical representation of input impedance versus frequency for a 
10 cage antenna optimized for VSWR < 2.0; 

FIG. 6A is a graphical representation of VSWR versus frequency for a cage 
antenna optimized for VSWR < 2.5 ; 

FIG. 6B is a graphical representation of directivity versus frequency for a cage 
antenna optimized for VSWR < 2.5; 
15 FIG. 7A is a perspective view of an exemplary cage monopole antenna 

embodiment in accordance with the presently disclosed technology, having dimensions a 
= 0.814mm, d= 2.2cm, w = 3.256mm, hj = 1.2cm, h 2 = 16cm and Z 0 = 50Q; 

FIG. 7B is a graphical representation of VSWR versus frequency for the 
exemplary cage monopole antenna of FIG. 7 A; 
20 FIG. 7C is a graphical representation of the input impedance versus frequency for 

the exemplary cage monopole antenna of FIG. 7 A; 

FIG. 7D is a graphical representation of the directivity versus frequency for the 
exemplary cage monopole antenna of FIG. 7 A; 

FIG. 8A is a perspective view of an exemplary sleeve-cage monopole antenna 
25 embodiment in accordance with the presently disclosed technology, having dimensions a 
= 0.814mm, d= 2.2cm, w = 3.256mm, h } = 1.2cm, h 2 = 16cm, r = 2.5cm, h'= 4cm, and Z 0 
= 50Q; 

FIG. 8B is a graphical representation of VSWR versus frequency for the 
exemplary sleeve-cage monopole antenna of FIG. 8 A; 
30 FIG. 8C is a graphical representation of the input impedance versus frequency for 

the exemplary sleeve-cage monopole antenna of FIG. 8 A; 
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FIG. 8D is a graphical representation of the directivity versus frequency for the 
exemplary sleeve-cage monopole antenna of FIG. 8A; 

FIG. 8E is a top view of the exemplary sleeve-cage monopole antenna of FIG. 

8A; 

5 FIG. 9A is a perspective view of an exemplary quadrifilar helical antenna 

embodiment in accordance with the presently disclosed technology, having dimensions a 
~ = 0.814mm, d = 2cm, w = 3.256mm, h } = 0.91cm, h 2 = 8.85cm, Z 0 = 50Q; 

FIG. 9B is a graphical representation of measured and computed VSWR versus 
frequency for the exemplary quadrifilar helical antenna of FIG. 9 A; 
10 FIG. 9C is a graphical representation of the measured and computed input 

impedance versus frequency for the exemplary quadrifilar helical antenna of FIG. 9 A; 

FIG. 9D is a graphical representation of the computed directivity (^=0) versus 
frequency for the exemplary quadrifilar helical antenna of FIG. 9 A; 

FIG. 1 OA is a perspective view of an exemplary sleeve helical antenna 
15 embodiment in accordance with the presently disclosed technology, having dimensions a 
= 0.814mm, d= 2cm, w = 3.256mm, hi = 0.91cm, h 2 = 8.85cm, Z 0 = 50Q; 

FIG. 1 0B is a graphical representation of measured and computed VSWR versus 
frequency for the exemplary sleeve helical antenna of FIG. 10A; 

FIG. 10C is a graphical representation of the measured and computed input 
20 impedance versus frequency for the exemplary sleeve helical antenna of FIG. 10A; 

FIG. 10D is a graphical representation of the computed directivity (^=0) versus 
frequency for the exemplary sleeve helical antenna of FIG. 10A; 

FIG. 1 1 A illustrates a curved wire helix for use in exemplary antenna technology 
of the present subject matter; 
25 FIG. 1 IB illustrates a curved wire loop for use in exemplary antenna technology 

of the present subj ect matter; 

FIG. 1 1 A illustrates a curved wire meander line for use in exemplary antenna 
technology of the present subject matter; 

FIG. 12 illustrates an arbitrary curved wire of radius a with source point, 
30 observation point, and unit vectors; 
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FIG. 13 illustrates an exemplary straight line segment approximation of a curved- 
wire axis; 

FIG. 14 illustrates exemplary geometric parameters for adjacent straight line 
segments; 

FIG. 15 represents a triangular basis function A n ; 

FIG. 16 illustrates an exemplary piecewise linear expansion of the current along a 
meander line; 

FIG. 17 represents a testing pulse n n ; 

FIG. 18 illustrates exemplary source and observation points on adjacent straight- 
wire segments; 

FIG. 19 illustrates exemplary source and observation points on different linear 
wire segments; 

FIG. 20 represents a composite triangle basis function with five constituent 
triangles; 

FIG. 21 represents an exemplary constituent triangle A* ; 
FIG. 22 represents an exemplary composite testing pulse TI p ; 

FIG. 23 represents an exemplary constituent testing pulse n£; 

FIG. 24 illustrates current expansion with composite basis functions; 
FIG. 25 A is a graphical representation of the real part of current versus degrees on 
wire loop antennas; 

FIG. 25B is a graphical representation of the imaginary part of current versus 
degrees on wire loop antennas; 

FIG. 26A is a graphical representation of the real part of current on the wire loop 
antenna of FIGS. 25 A and 25B without composite basis function at the source; 

FIG. 26B is a graphical representation of the imaginary part of current on the wire 
loop antenna of FIGS. 25 A and 25B without composite basis function at the source; 

FIG. 27 is a graphical representation of mode 2 current on one arm of four arm 
Archimedian spiral antenna having spiral constant 0.02 A,, arm length \2X and wire radius 
0.006A.; 
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FIG. 28 is a graphical representation of the magnitude of mode 2 current on one 
arm of four arm Archimedian spiral antenna having spiral constant 0.021, arm length 121 
and wire radius 0.0061; 

FIG. 29 is a graphical representation of the magnitude of mode 2 current on one 
5 arm of four arm Archimedian spiral antenna having spiral constant 0.02 A,, arm length 12 A, 
and wire radius 0.0061 with data from [2] Champagne et al. for comparison; 

FIG. 30 A is a perspective view of an exemplary helical geometry; 

FIG. 3 0B represents additional exemplary definitions of helical parameters; 

FIG. 31 is a graphical representation of current versus arc displacement on a helix 
10 illuminated by a plane wave, E = (jccostf + zsm0)e- Mxs ^ zcos0) , 

(0= 45°, L = 0.51, v = 10, a= 20°, a = 0.00051); 

FIG. 32 is a graphical representation of current versus arc displacement on a helix 
illuminated by a plane wave, E = (xcos0 + zsin0)e- Jk{xsinO ' zcosO) , 

(0= 45°, L = 0.351, v = 10, a= 20°, a = 0.00051); 
15 FIG. 33 is a graphical representation of current versus arc displacement on a helix 

illuminated by a plane wave, E = ze~ jkx , (L = 21, v = 50, a- 20°, a = 0.0005?i); 

FIG. 34 is a graphical representation of current versus arc displacement on a helix 
driven by delta-gap source (L = 0.25 A, v = 25, <z= 20°, a = 0.0005A,); 

FIG. 35 is a graphical representation of current versus arc displacement on a helix 
20 driven by delta-gap source (L = X, v = 25, a= 20°, a = 0.00051); 

FIG. 36 is a graphical representation of current versus arc displacement on a helix 
driven by delta-gap source (L = 0.51, Ch = 0.11, a= 12.5°, a = 0.0051); 

FIG. 37 is a graphical representation of VSWR versus frequency for different 
pitch angles for an antenna helix; 
25 FIG. 38 is a graphical representation of VSWR versus frequency for different 

wire radius values for an antenna helix; 

FIG. 39 provides a block diagram of exemplary steps in a process for efficient 
evaluation of antennas; 

FIG. 40A is a graphical representation of VSWR versus frequency for the 
30 exemplary antenna of FIG. 40B; 
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FIG. 40B is a perspective view of an exemplary straight wire antenna with two 
parasites; 

FIG. 41 A is a graphical representation of VSWR versus frequency for the 
exemplary antenna of FIG. 4 IB; 
5 FIG. 41 B is a perspective view of an exemplary straight wire antenna with four 

parasites; 

FIG. 42 A is a graphical representation of VSWR versus frequency for the 
exemplary antenna of FIG. 42B; 

FIG. 42B is a perspective view of an exemplary helix antenna with two straight 
10 wire parasites; 

FIG. 43 A is a graphical representation of VSWR versus frequency for the 
exemplary antenna of FIG. 43B; 

FIG. 43 B is a perspective view of an exemplary helix antenna with four straight 
wire parasites; 

15 FIG. 43 C is a graphical representation of directivity versus frequency for the 

exemplary antenna of FIG. 43B; 

FIG. 43 D is a graphical representation of directivity in the H-plane versus § for 
the exemplary antenna of FIG. 43B; 

FIG. 44 A is a graphical representation of VSWR versus frequency for the 
20 exemplary antenna of FIG. 44B; 

FIG. 44B is a perspective view of an exemplary helix antenna with two helical 
parasites; 

FIG. 44C is a graphical representation of input impedance versus frequency for 
the exemplary antenna of FIG. 44B; 
25 FIG. 44D is a graphical representation of directivity versus frequency for the 

exemplary antenna of FIG. 44B; 

FIG. 45 A is a graphical representation of VSWR versus frequency for the 
exemplary antenna of FIG. 45B; 

FIG. 45B is a perspective view of an exemplary helix antenna with inner and 
30 outer helical parasites; 
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FIG. 46A is a perspective view of an exemplary antenna base portion and 
representative cylinder around which the coils of a triple helix antenna are wound; 

FIG. 46B is a perspective view of an exemplary triple helix antenna; 

FIG. 46C is a graphical representation of VSWR versus frequency for a triple 
5 helix antenna (such as one depicted in FIGS. 46 A and 46B) compared with a single helix 
antenna; 

FIG. 47A is a graphical representation of VSWR versus frequency for the 
exemplary antenna of FIG. 47B compared with a single helix and a triple helix antenna; 

FIG. 47B is a perspective view of an exemplary triple helix antenna with four 
10 straight wire parasites; 

FIG. 47C is a graphical representation of directivity versus frequency for a triple 
helix antenna; 

FIG. 47D is a graphical representation of directivity versus frequency for a triple 
helix antenna with parasites, such as the one illustrated in FIG. 47B; 
15 FIG. 48 A is a graphical representation of VSWR versus frequency for cage 

monopole antenna optimization for VSWR < 2.5; 

FIG. 48B is a graphical representation of directivity versus frequency for cage 
monopole antenna optimization for VSWR < 2.5; 

FIG. 49A is a graphical representation of VSWR versus frequency for cage 
20 monopole antenna optimization for VSWR < 2.0; and 

FIG. 49B is a graphical representation of directivity versus frequency for cage 
monopole antenna optimization for VSWR < 2.0. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
25 DESIGN PROCEDURE: 

We modeled and measured the properties of a so-called cage monopole. The cage 
monopole shown in FIG. 1 A consists of four vertical straight wires connected in 
parallel and driven from a common stalk at the ground plane. The ground plane in this 
model is assumed to be of infinite extent to facilitate analysis. The parallel straight wires 
30 are joined by crosses constructed of brass strips. The strip width was selected to be 
electrically equivalent to the wire radius for the purpose of modeling the structure. 
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Compared to a single, thin, straight wire, the cage structure with multiple wires has a 
lower peak voltage standing wave ratio (VSWR) over the band. This is important since a 
structure which has a comparatively small VSWR over a band is more amenable to 
improvements in bandwidth with the addition of other components such as loads or 
5 parasites than is the common single-wire monopole with higher VSWR. 

Next we add four parasitic straight wires of equal height (h) and distance (r) from 
the center of the cage to create the so-called "sleeve-cage monopole" of FIG. 2 A. The 
genetic algorithm of (D.L. Carroll, "A FORTRAN Genetic Algorithm Driver", Univ. of 
Illinois, Urbana, IL, http://www.stafF.uiuc.edxi/-carroll/ga.html) is used to determine the 

10 optimum distance and height of these parasitic straight wires. In this example the fitness 
value assigned to each antenna in the optimization process is the bandwidth ratio defined 
by f 2 l f x , where fi and fx are, respectively, the largest and smallest frequencies 
between which the VSWR is 3.5 or less. We interpolate the impedance matrix with 
respect to frequency in order to increase the speed of the optimization process. 

15 To design antennas that are smaller, we turn our attention to the normal mode helix, 

since, for operation about a given frequency, it can be made shorter than the vertical whip 
by adjustment of the pitch angle. Also, we observe a decrease in the peak VSWR when 
additional filaments are added to the helix driven from a central straight wire. Generally, 
a normal mode helix will exhibit electrical properties similar to those of a straight wire 

20 having the same wire length, though the peak VSWR for the helix is usually greater. The 
quadrifilar helix of FIG. 3 A whose height is 9.8 cm can be used in the same bands as 
a cage monopole of height about 14 cm. Thus, the total height of the antenna can be 
reduced by 30% with the 42° pitch angle. Parasitic straight wires of optimum height and 
distance are added to create what we call the "sleeve helical monopole" shown in FIG. 

25 4A. 

RESULTS AND DISCUSSION: 

As one can see from the VSWR data, good agreement is achieved between 
predictions computed by means of our numerical techniques and results measured on a 
30 model mounted over a large ground plane. The frequency range over which data are 

presented is dictated by the frequencies over which our ground plane is electrically large. 
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The slight discrepancies in the computed and measured results are attributed to 
imprecision in the construction of the antennas. The predicted results of bandwidth and 
VSWR of each antenna are summarized in Table 1 below. 



Structure 


\/o\i/p 
VjWK 


D W IVdllU 


D W /0 


Prpni ipnp\/ Rutictp 
riCUUCllL>Y IVCUlgC - 


(cm) 


Width 
(cm) 




<5.0 


1 1.7:1 


312 


300-3500 


17.2 


2.2 


<3.5 


3:1 


115 


950-2850 


Sleeve-cage 
monopole 


<5.0 


5.2:1 


185 


315-1650 


17.2 


5 


<3.5 


4.4:1 


163 


350-1550 


Quadrifilar helix 


<5.0 


5.8:1 


199 


475-2750 


9.8 


2 


<3.5 


1.6:1 


47 


500-800 


Sleeve helix 


<5.0 


3.9:1 


147 


475-1850 


9.8 


6 


<3.5 


3.5:1 


134 


500-1750 



Table 1 



5 We point out that when the parasitic elements are added to each structure, the 

bandwidth ratio increases for the VSWR < 3.5 requirement. However, outside of this 
frequency range the VSWR is worse than that of the antenna without parasites. In other 
words, VSWR has, indeed, been improved markedly over the design range but at a 
sacrifice in performance outside the range, where presumably the antenna would not be 

10 operated. Also, notice that the deep nulls in the directivity at the horizon for the cage and 
the quadrifilar helix structures have been eliminated with the addition of the parasites. 
Thus the directivity is improved in the band where on the basis of VSWR this antenna is 
deemed operable, although there was no constraint on directivity specified in the 
objective function. 

15 CAGE ANTENNAS OPTIMIZED FOR BANDWIDTH 

Design Method: The cage antenna is depicted in FIG. 7A where one sees four 
vertical wires joined to the feed and stabilized by thin brass strips of width w. The strips 
are treated as wires of radius a= w/4 . The GA of (CARROLL, D.L.: 'Chemical Laser 
Modeling with Genetic Algorithms', AIAA Journal Feb. 1996, 34, (2), pp. 338-346) is 

20 applied to optimize the diameter (d) of the cage structure and the length Q12) of the 

wires in the cage. Each function evaluation consists of numerically solving the electric 
field integral equation for the cage geometry (having dimensions chosen by the GA) over 
the band of interest. Candidate antennas are given a fitness score equal to the bandwidth 
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ratio f h /f x , where fi is the lowest and fo is the highest frequency of operation over a 

band where the VSWR meets the design goal. 

Results: The antenna whose characteristics are represented in FIGS. 5A and 5B is 
optimized for a design goal of VSWR < 2.0 over the frequency band 500 to 1600 MHz. 

5 The GA picks the parameter d from a range of 1 cm to 5 cm with a resolution of 0. 1 3 cm 
(5 bits, 32 possibilities). The range specified for parameter h 2 is 8 cm to 12 cm with a 
resolution of 0.27 cm (4 bits, 16 possibilities). The GA converges to an optimum solution 
after three generations with five antennas per generation. A sensitivity analysis reveals 
that antenna input characteristics change only modestly with small geometric variation. 

10 The directivity of this cage antenna for <f> = 0° and <f> = 75°,90° is above 4 dBi over the 
entire band. The properties of this antenna and those of Nakano's helical monopole 
(NAKANO, H., IKEDA, N., WU, Y., SUZUKI, R., MIMAKI, H., and YAMAUCHI, J.: 
'Realization of dual-frequency and wide-band VSWR performances using normal-mode 
helical and inverted-F antennas', IEEE Trans. Antennas Propagate June 1998, 46, (6), 

15 pp. 788-793), which is designed to operate with VSWR < 2.0 in two frequency bands, are 
listed in Table 2 for comparison. 



Structure 


Cage 
(Figure 6) 


HX-MP 


Cage 
(Figure 7) 


Sleeve Dipole 


VSWR 


<2 


<3.5 


<2.5 


<2.5 


Bandwidth Ratio 


2.6 


1.7 


5.4 


1.8 


f (MHz) 


575-1500 


627-1048 


210-1130 


225-400 


Height (cm) 


10.3 


19.8 


26.55 


51 


Width (cm) 


4.7 


0.47 


8.2 


13 


Wire radius (mm) 


0.814 


0.3 


3.175 


14.3 



Table 2 



The antenna whose characteristics are represented in FIGS. 6 A and 6B is 
20 optimized for a design goal of VSWR < 2.5 in the frequency range 200 to 1200 MHz. 
This range is chosen for comparison of the cage antenna to the open sleeve dipole of 
(KING, H.E. 5 and WONG, J.L.: 'An experimental study of a balun-fed open-sleeve 
dipole in front of a metallic reflector', IEEE Trans. Antennas Propagate March 1972, 20, 
(2), pp. 201-204) which operates over the frequency range 225 to 400 MHz. The GA is 
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allowed to choose parameter d from 1 cm to 10 cm with a resolution of 0.6 cm (4 bits, 16 
possibilities). The parameter Pi2 is selected from 20 cm to 25 cm with a resolution of 0.33 
cm (4 bits, 16 possibilities). An optimum result is reached after 1 1 generations with five 
antennas per generation. This cage monopole is not useful over the entire frequency range 
5 for which its VSWR is less than 2.5 since there is a null in the directivity within this 
range. It is operable over a 3.6:1 bandwidth for VSWR less than 2.5 and directivity 
greater than 0 dBi. In Table 2 are listed the properties of the cage antenna together with 
those of the sleeve dipole. 

CAGE MONOPOLE AND SLEEVE-CAGE MONOPOLE 

10 The cage monopole shown in FIG. 7 A consists of four vertical straight wires 

connected in parallel and driven from a common wire which is the extension of the center 
conductor of a coaxial cable protruding from the ground plane. The ground plane in this 
model is assumed to be of infinite extent in the analysis of the structure. The parallel 
straight wires are joined by crosses constructed of brass strips. The strip width w was 

15 selected to be electrically equivalent to the wire radius a for the purpose of modeling the 
structure (w=4 a ) (CM. Butler, "The equivalent radius of a narrow conducting strip," 
IEEE Trans. Antennas Propagate vol. AP-30, pp.755-758, July 1982). Compared to a 
single, thin, straight wire, the cage structure with multiple wires has a lower peak VSWR 
over the band as seen in FIG. 7B. This is important since a structure which has a 

20 comparatively small VSWR over a band is more amenable to improvements in bandwidth 
with the addition of other components such as loads or parasites than is the common 
single-wire monopole with higher VSWR. 

Four parasitic straight wires of equal height (h) and radial distance (r) from the 
center line of the cage are added to create the so-called "sleeve-cage monopole" of FIG. 

25 8A. The genetic algorithm of (D.L. Carroll, "Chemical Laser Modeling with Genetic 
Algorithms ,'M,kL4 Journal, vol. 34, no. 2), pp. 338-346, Feb. 1996) is used to determine 
optimum values of h and r for given design goals. An objective function evaluation for 
one antenna in the GA population involves numerically solving the electric field integral 
equation for many frequencies within the band of interest. Since this must be done for 

30 many candidate antennas, it is advantageous to interpolate the integral equation 

impedance matrix elements with respect to frequency (E.H. Newman, "Generation of 
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wide-band data from the method of moments by interpolating the impedance matrix," 
IEEE Trans. Antennas Propagate vol. AP-36, pp. 1820-1824, Dec. 1988). Each 
candidate structure is assigned a fitness value based on its electrical properties. A simple 
fitness value used here is the antenna bandwidth ratio which measures the performance of 
5 the antenna over a frequency band of interest denoted by [f A , f B ] . The bandwidth ratio 
for a particular antenna is considered a function of its geometry and is computed from 

F(Kr) =A 
J\ 

where f x = min(/) such that VSWR(f) < limit 

MfaJ B ] 

and f 2 = max(/) such that VSWR(f) < limit for all / e [/J , f 2 ] . 

MAJb] 

10 Another viable fitness value is the percent bandwidth defined here as 

%BW = 1Q0^J±. 

QUADRIFILAR HELIX AND SLEEVE HELIX 

To design low profile antennas, we turn our attention to the normal mode helix, 
15 since, for operation about a given frequency, it can be made shorter than the vertical whip 
by adjustment of the helix pitch angle. Generally, a normal mode helix will exhibit 
electrical properties similar to those of a straight wire having the same wire length, 
though the peak VSWR for the helix is usually greater. The helix exhibits vertical 
polarization as long as it operates in the normal mode. There is a decrease in the peak 
20 VSWR, relative to that of a single-wire helix, when additional helical filaments are added 
to one driven from a central straight wire. The quadrifilar helix of FIG. 9A whose 
height is 9.8 cm can be used in the same bands as a cage monopole of height about 14 
cm. Thus, the total height of the antenna can be reduced by 30% with the 42° pitch 
angle. Parasitic straight wires of optimum height and distance are added to create what- 
25 we call the "sleeve helical monopole" shown in FIG. 10A. Most integral equation 

solution techniques for the helix are, in general, more computationally expensive since 
these require many basis functions to represent the vector direction of the current along 
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the meandering wire. A solution procedure which uncouples the representation of the 
geometry from the representation of the unknown current is developed in (S.D. Rogers 
and CM. Butler, "An efficient curved-wire integral equation solution technique," 
submitted to IEEE Trans. Antennas Propagat.) and is used here to reduce the time in 
5 optimization of antennas with curved wires. 
RESULTS 

As one can see from the VSWR data, good agreement is achieved between 
predictions computed by means of numerical techniques (S.D. Rogers and CM. Butler, 
"An efficient curved-wire integral equation solution technique," submitted to IEEE 

10 Trans. Antennas Propagat) and results measured on a model mounted over a large 

ground plane. The frequency range over which our experiments are conducted is dictated 
by the frequencies over which the ground plane is electrically large. Of course, the 
dimensions of the antenna may be scaled for use in other bands. The slight discrepancies 
in the computed and measured results are attributed to the difficulty in building the 

15 antenna to precise dimensions. However, a sensitivity analysis reveals that the antenna 
performance changes minimally with small variations in geometry. The reflection 
coefficient is measured at the input of the coaxial cable driving the monopoles and of a 
shorted section of coaxial line having the same length. Applying basic transmission line 
theory to these data, one can determine the measured input impedance of the antenna 

20 with the reference "at the ground plane." All VSWR data is for a 50Q system. As the 

feed point properties of the various antennas are evaluated, we must also keep in mind the 
radiation properties of the antenna, so computed directivity is included herein. The 
predicted results of bandwidth and VSWR of each antenna are summarized in Table 1. 
We point out that, when the parasitic elements are added to each structure, the 

25 bandwidth ratio increases for the VSWR < 3.5 requirement. However, outside of this 
frequency range the VSWR is worse than that of the antenna without parasites. In other 
words, VSWR has, indeed, been improved markedly over the design range but at a 
sacrifice in performance outside the range, where presumably the antenna would not be 
operated. Also, notice that the deep nulls in the directivity at the horizon for the cage and 

30 the quadrifilar helix structures have been eliminated with the addition of the parasites. 
Thus the directivity is improved in the band where, on the basis of VSWR, this antenna is 



deemed operable, although there was no constraint on directivity specified in the 
objective function. 

A summary and comparison of the results for the various antenna structures 
represented in FIGS. 7 A, 8 A, 9 A and 10A as well as Nakano's Helix Monopole and a 
SDNfCGARS dipole antenna that was developed and produced by ITT for the Army is 



included in the following Table 3. 



Structure 


VSWR 


BW Ratio 

A 
f 2 


BW% 

100 (/,"/ 2 ) 
V/1/2 


Frequency Range 
(MHZ) 


Height 
(cm) 


Width 
(cm) 


Cage monopole 


<3.5 


3:1 


115 


950-2850 


17.2 


2.2 


Sleeve-cage 
monopole 


<3.5 


4.4:1 


163 


350-1550 


17.2 


5 


Quadrifilar 
helix 


<3.5 


1.6:1 


47 


500-800 


9.8 


2 


Sleeve helix 


<3.5 


3.5:1 


134 


500-1750 


9.8 


6 


SINCGARS 
Antenna 


<3.5 


2.9 


112 


30-88 


280 


2 


Nakano's Helix 
Monopole 


<3.5 


1.7 


52 


627-1048 


19.8 


0.4 



Table 3 



Additional results are now presented for the antenna of FIG. 7A with optimization 
for VSWR of less than 2.5. An antenna is tested having the following parameters: a = 
3.175mm, d= 7.6cm 5 w = 1.27 cm, hj = 2.55cm and h 2 = 22.95cm. As seen in the graph 
of FIG. 48A, a VSWR of less than 2.5 is achieved over a frequency band of 212-1 155 
MHz, resulting in a bandwidth ratio of 5.5: 1. FIG. 48B shows the directivity versus 
frequency for the same antenna for different angles of theta (0). A summary and 
comparison of results for the optimization with VSWR<2.5 described above is listed in 



Table 4 below. 



Structure 


VSWR 


BW Ratio 
fx 

f 


BW% 

100 uwp 


Frequency Range 
(MHZ) 


Height 
(cm) 


Width 
(cm) 


King's open 
sleeve dipole 


<2.5 


1.77 


58.3 


225-400 . 


51 


13 


NTDR Antenna 


<2.5 


2.0 


70 


225-450 


200 


6.4 


Cage monopole 


<2.5 


3.7 


139 


212-775 


25.5 


7.6 



Table 4 



Additional results are now presented for the antenna of FIG. 7A with optimization 
for VSWR of less than 2.0. An antenna is tested having the following parameters: a = 
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0.814mm, d= 4.8cm, w = 3.256cm, hi = 1cm and h 2 = 9cm. Results are presented in 
FIG. 49A for the described cage monopole antenna as well as for thin straight wire 
antenna and a fat straight wire antenna. As seen in the graph of FIG. 49 A, a VSWR of 
less than 2.0 is achieved over a frequency band of 575-1 500 MHz, resulting in a 



5 bandwidth ratio of 2.6: 1 . FIG. 49B shows the directivity versus frequency for the same 
antenna for different angles of theta (0). A summary and comparison of results for the 
optimization with VSWR<2.0 described above is listed in Table 5 below. 



Structure 


VSWR 


BW Ratio 

A 
h 


BW% 

100 c/;-/2) 

V/./2 


Frequency 
Range (MHZ) 


Height 
(cm) 


Width 
(cm) 


a 

(mm) 


Nakano's Helix 
Monopole 


<2.0 


1.14 


13.44 


662-757 


19.8 


0.4 


0.015, 
0.003 


<2.0 


1.05 


5.78 


957-1014 


Cage monopole 


<2.0 


3.7 


139 


212-775 


25.5 


7.6 


0.814 



Table 5 



Conclusions from the above numerical results include recognition that cage 
10 structures can be optimizes for lower VSWR, parasites of optimum size and placement 
improve VSWR of driven antenna, helical elements reduce height at sacrifice of 
bandwidth, and wire radius is an important parameter. 



The following is a detailed description (including documentary references, a list 
15 for which is provided after the detailed description) of an exemplary efficient curved- wire 
integral equation solution technique as may be practiced in accordance with the subject 
invention. 

AN EFFICIENT CURVED-WIRE INTEGRAL EQUATION SOLUTION 
20 TECHNIQUE: 

Computation of currents on curved wires by integral equation methods is often 
inefficient when the structure is tortuous but the length of wire is not large relative to 
wavelength at the frequency of operation. The number of terms needed in an accurate 
piecewise straight model of a highly curved wire can be large yet, if the total length of 
25 wire is small relative to wavelength, the current can be accurately represented by a simple 
linear function. In embodiments of the present invention, a new solution method for the 
curved- wire integral equation is introduced. It is amenable to uncoupling of the number 
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of segments required to accurately model the wire structure from the number of basis 
functions needed to represent the current. This feature lends itself to high efficiency. The 
principles set forth can be used to improve the efficiency of most solution techniques 
applied to the curved-wire integral equation. New composite basis and testing functions 
5 are defined and constructed as linear combinations of other commonly used basis and 
testing functions. We show how the composite basis and testing functions can lead to a 
reduced-rank matrix which can be computed via a transformation of a system matrix 
created from traditional basis and testing functions. Supporting data demonstrate the 
accuracy of the technique and its effectiveness in decreasing matrix rank and solution 

10 time for curved- wire structures. 
I. INTRODUCTION: 

Numerical techniques for solving curved- wire integral equations [1] may involve 
large matrices, often due primarily to the resources needed to model the structure 
geometry rather than due to the number of basis functions needed to represent the 

15 unknown current. This is obviously true when a subdomain model is used to approximate 
a curvilinear structure in which the total wire length is small compared to the wavelength 
at the frequency of operation. Usually the number of segments needed in such a model is 
dictated by the structure curvature rather than by the number of weighted basis functions 
needed in the solution method to represent the unknown current. There is a demand for a 

20 general solution technique in which the number of unknowns needed to accurately 

represent the current is unrelated to the number of straight segments required to model 
(approximately) the meandering contour of the wire and the vector direction of the 
current. In recent years attention in the literature has been given to improving the 
numerical efficiency of integral equation methods for curved- wire structures [2]-[13]. For 

25 the most part, presently available techniques incorporate basis functions defined on 
circular or curved wire segments. The authors of [2] define basis and testing functions 
along piecewise quadratic wire segments and achieve good results with fewer unknowns 
than would be needed in a piecewise straight model of a wire loop and of an Archimedian 
spiral antenna. Others introduce solution techniques for structures comprising circular 

30 segments that numerically model the current specifically on circular loop antennas [3], 
[4], An analysis of general wire loops is presented in [5], where a Galerkin technique is 
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employed over a parametric representation of a superquadric curve. In [6] arcs of 
constant radii are employed to define the geometry of arbitrarily shaped antennas from 
which is developed a technique for analyzing helical antennas. Other methods which 
utilize curved segments for subdomain basis and testing functions are available [7]-[10]. 
5 There are several advantages inherent in techniques in which basis and testing 

functions are defined over curved wire segments. Geometry modeling error can be made 
small and solution efficiency can be increased since to "fit" some structural geometries 
fewer curved segments are needed than is feasible with straight segments. Although these 
techniques are successful, they suffer disadvantages as well. First, the integral equation 

10 solution technique must be formulated to account for the new curved-segment basis and 
testing functions. This means that computer codes must be written to take advantage of 
the numerical efficiency of these new formulations incorporating the curved elements. A 
second disadvantage of curvilinear basis function modeling is that they fit one class of 
curve very well but are not well suited to structures comprising wires of mixed curvature. 

15 That is, circles fit loops and helices well but not spirals. Clearly, when a given structure 
comprises several arcs of different curvatures, the efficiency of methods employing a 
single curved-segment representation suffers. Elements like the quadratic segment or 
the arc-of-constant-radius segment increase the complexity of modeling. The third 
disadvantage of these techniques is that, for many structures, they do not lead to complete 

20 uncoupling of the number of the unknown current coefficients from the number of 
segments needed to model the structure geometry. For example, several quadratic 
segments or arcs, with one weighted unknown defined on each, would be required to 
model the geometry of one turn of a multiturn helix, yet the current itself may be 
represented accurately in many cases by a simple linear function over several 

25 turns. 

In this description, an efficient method for solving for currents induced on curved- 
wire structures is presented. The solution method is based on modeling the curved wire 
by piecewise-straight segments but the underlying principles are general and can be 
exploited in conjunction with solution procedures which depend upon other geometry 
30 representations, including those that use arcs or curves. It is ideal for multi-curvature wire 
structures [12], [13]. The improved solution technique depends upon new basis and 
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testing functions which are defined over more than two contiguous straight- wire 
segments. Composite basis functions are created as sums of weighted piecewise linear 
functions on wire segments, and composite testing functions compatible with the new 
basis functions are developed. The new technique allows one to reduce the rank of 
5 the traditional impedance matrix. We show how the matrix elements for a reduced-rank 
matrix can be computed from the matrix elements associated with a traditional integral 
equation solution method. Of paramount importance is the fact that the number of 
elements employed to model the geometric features of the structure is unrelated to the 
number of unknowns needed to accurately represent the wire current. 

10 The concept of creating a new basis function as a linear combination of other 

basis functions is used in [14] for a multilevel iterative solution procedure for integral 
equations. Perhaps the composite basis function defined herein can be thought of as a 
"coarse level" basis function in multilevel terminology, although the method described in 
this specification is not related to the so-called multilevel or multigrid theory of [14]-[16]. 

15 The improved solution technique requires fewer unknowns than the traditional 

solution to represent the current on an Archimedian spiral antenna. Results comparable to 
those presented in [2] are achieved for the spiral. The improved technique also allows one 
to significantly reduce the number of unknowns required to solve for the current on wire 
helices. Specifically, the results of a convergence test show that the current on a helix can 

20 be modeled accurately with the same number of unknowns needed for a "similar" straight 
wire even though the helix has a large number of turns. 
II. INTEGRAL EQUATION FOR GENERAL CURVED WIRES: 

In this section we present the integro-differential equation governing the electric 
current on a general three-dimensional curved or bent wire. Examples are the wire loop, 

25 the helix, and the meander line shown respectively in FIGS. 1 1 A, 1 IB and 1 1C. The wire 
is assumed to be a perfect electrical conductor and to be thin which means that the radius 
is much smaller than the wavelength and the length of wire. Under these thin-wire 
conditions the current is taken to be axially directed, circumferentially invariant, and zero 
at free ends. The equation governing the total axial current I(s)s on the thin curved wire 

30 is 
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-jJLh 1 \l(s'#'-tiC(s 9 s'W+— {—I(s')K(s 9 s')ds\ = -E i (;s)-* 9 seC (1) 
4xk [ * ds Ids' J 

in which C is the wire axis contour, s denotes the arc displacement along C from a 
reference to a point on the wire axis, and s is the unit vector tangent to C at this point. 
The positive sense of this vector is in the direction of increasing s . K{s 9 s') is the kernel 
5 or Green's function, 

K{sJ)~ f—df 9 (2) 

27t_{ R 

in which R is the distance between the source and observation points on the wire surface, 
andE' (s) is the incident electric field which illuminates the wire, evaluated in (1) on the 
wire surface at arc displacement s . Geometric parameters for an arbitrary curved wire are 

10 depicted in FIG. 12. 

III. TRADITIONAL SOLUTION TECHNIQUE: 

The new solution method proposed in this specification can be viewed as an 
improvement to present methods. In fact, employing the ideas set forth in Section IV, one 
can modify an existing subdomain solution method to render it more efficient for solving 

15 the curved-wire integral equation. Hence, the new method is explained in this 

specification as an enhancement of a method that has proved useful for a number of 
years. The method selected for this purpose is based on modeling the curved wire as an 
ensemble of straight-wire segments, with the unknown current represented as a linear 
combination of triangle basis functions and testing done with pulses. In this section this 

20 method is outlined as a basis for the explanation of the new method in Section IV. 

The first step in modeling a curved wire is to select points on the wire axis and 
define vectors r 0 ,r, ,...,r from a reference origin to the selected points. The curved wire 

is modeled approximately as an ensemble of contiguous straight- wire segments joining 
these points (cf. FIG. 13). The arc displacement along the axis of the piecewise linear 
25 approximation of C is measured from the reference point labeled r 0 . The arc 

displacement between r 0 and the n th point located by r„ is /„ . A general point on the 
piecewise-straight approximation of the wire axis is located alternatively by means of the 
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vector r and by the arc displacement / from the reference to the point. Various 
geometrical parameters describing the wire can be expressed in terms of the 
vectors locating the points on the wire axis. The unit vectors along the directions of the 
segments adjacent to the point r p shown in FIG. 14 are given by 

5 U-^V^ 0) 

v 

i = r p +i ~ r p .(4) 

where 

10 A p+ =1^,-^1- (6) 

The midpoint of the straight-wire segment joining r p and r p±1 is located by 



I 

t p±- 2 



r +r . (7) 



In order to emphasize the fact that the model is now a straight wire segmentation 
of the original curved wire, s in (1) is replaced by /, the arc displacement along the axis of 
15 the straight wire model. With this notation and subject to the piecewise straight wire 
approximation, Eq. (1) becomes 

-j-rj\k 2 f/(/')rk(/,ryr + 4 &mK(i,r)di\ = -e'(/)-u e l (8) 

4^fe[ l dl[dV J 

where L is the piecewise straight approximation to C . 

In a numerical solution of the integral equation for a curved wire structure, the 
20 (vector) current is expanded in a linear combination of weighted basis functions defined 
along the straight- wire segments. Even though they can be any of a number of functions, 
those employed here, for the purpose of illustration in this specification, are chosen to be 
triangle functions with support over two adjacent segments. Thus the current may be 
approximated by 
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A n (/)i n (0 (9) 

n=\ 
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in which the triangle function A„ about the n' h point on the segmented wire, as depicted in 
FIG. 15, is defined by 

'/-/ 



A (/) = 



n-\ 



n—l n 



I -I 

n+l 



(10) 



WSJ 



where the unit vector I is defined in terms of the unit vectors associated with the 



th 

5 segments adjacent to the n point: 



1 = 



(11) 



N is the number of basis functions and unknown current coefficients I„ in the finite series 
approximation (9) of the current. TV unknowns are employed to represent the current on a 

10 wire having two free endpoints and modeled by N +1 straight- wire segments. In this 

traditional solution technique described here, TV must be large enough to accurately model 
the geometric structure and vector direction of the current, even if a large number of 
unknowns is not required to approximate the current 1(1) to the accuracy desired The 
triangle basis functions overlap as suggested in FIG. 16 so an approximation with N 

15 terms incorporates, at most, +1 vector directions of current on the wire. These point- 
by-point directions of current on a curved wire must be accounted for accurately by the N 
+1 unit vectors, yet TV piecewise linear basis functions may be far more than may be 
needed to accurately represent the current 1(1). 

Testing the integro-differential equation is accomplished by taking the inner 

20 product of (8) with the testing function 



n (/) = 



! » /e < 7 . > l J 

m— m+ 

0, otherwise 



(12) 



28 



10 



depicted in FIG. 17 for m '- 1, 2,. . N. The inner product of this testing pulse with a 
function of the variable / is defined by 



f,n) = \~f(!)di. 



(13) 

Expanding the unknown current / with (9) and taking the inner product of (8) with (12) 
for m -\, 2,..., N yield a system of equations written in matrix form as 

[ZJU>[TJ (14) 

where 



mn n m 



z=-J 



(A 1 1 + A 1 1 ) f " A (T)K(R )dV 
+ (A 1 1 + A 1 J )f"'A (l')K(R )dV 

n 



A 



— ("+'K(R )dV> 
A J - 

is an element of the Nx N impedance matrix with 

" r 



4a 2 sin 


2 ^+(i -i) 2 

2 


r -r' 

m 


2 + a\ 



I and / ' on same segment 

m 

otherwise 



and 



R = 

m± 



4a 2 sin 2 -^ + (/ -O 2 , 



(15) 



(16) 



r -r 

m± 



,2 + a 2 , 



I and / ' on same segment 

m± 



otherwise 



(17) 



When the source (r' or /') and observation (r or / = (/ ,1 J) points reside on the same 
1 5 straight wire segment of radius a, as in FIG. 1 8 the exact kernel given by 
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,e- JkR 



"w-khr* (18) 



is used. Otherwise for source and observation points on different straight- wire segments 
(cf. FIG. 19), the exact kernel is approximated by the so-called reduced kernel, 



, -jkR 



K(i,n=-—. (19) 

K 

The approximation below, which is excellent when the segment lengths are small 
compared with the wavelength, is employed in arriving at the first two terms of (15): 



(i ( /).f(/) 5 n m( /))*f(/ m ). 



1 * l 

—A 1 + -A 1 

2 m ~ m ~ 2 



m+ m+ 



(20) 



The same approximation can be used to compute the elements of the excitation column 
vector, 



10 F = -£'(/)• i,n >-£'(/ ) 



1 « 1 

—A 1 +-A 1 

2 m- m- 2 m+ m+ 



(21) 



where E'(7^) is the known incident electric field at point l m on the wire. Of course, if 
desired the left hand side of (21) can be evaluated numerically in those situations in 
which the incident field varies appreciably over a subdomain. We also point out that 
testing with pulses allows one to integrate directly the second term on the left side of (8). 

15 The derivative of the piecewise linear current in (8) leads to a pulse doublet (for charge) 
over two adjacent straight wire segments. These operations on the second term in the left 
side of (8) lead to the last four integrals in (15). 
IV. IMPROVED SOLUTION TECHNIQUE 

In this section a new technique for solving the curved- wire integral equation is 

20 presented. It is very efficient for tortuous wires on which the actual variation of the 

current is modest, a situation which often occurs when the length of wire in a given curve 
is small relative to wavelength, regardless of the degree of curvature. Composite basis 
and testing functions are introduced as an extension of the functions of the traditional 
solution method outlined in Section III. The composite basis function serves to uncouple 

25 the number of straight segments needed to model the curved- wire geometry and the 
vector direction of the current from the number of unknowns needed to accurately 
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represent the current on the wire. This new basis function is a linear combination of 
appropriately weighted generic basis functions, e.g., basis functions (9) in the traditional 
method outlined in Section III, and is defined over a number of contiguous straight 
segments. This new basis function is referred to as a composite basis function since it is 

5 constructed from others. Even though the solution method can incorporate any number of 
different generic basis and testing functions, the piecewise linear or triangle basis 
function and the pulse testing function are adopted here to facilitate explanation. Also, 
this pair leads to a very efficient and practicable solution scheme. 

The notion of a composite triangle made up of constituent triangles is suggested 

10 in FIG. 20. For simplicity in illustration, the composite triangle is shown over a straight 
line though in practice it would be over a polygonal line comprising straight-line 
segments, which approximate the curved wire axis. The q th composite vector triangle 
function can be constructed as 



MO* =2>;a;(/)i; (22) 

1 5 in which A * is the i' h constituent triangle defined by 

l-l" 



A ? (/)l" = < 



A" 

i- 



/*-/ 

— 1', leWJ 



(23) 



A" 

i+ 



and illustrated in FIG. 21 . When q is used as a superscript it identifies a parameter related 
to the q th composite triangle function. The constitutive elements of the q th composite basis 

function are denoted by the subscript /. The parameter h q is the weight or magnitude of 



20 the / constituent triangle within A^ . These weights are functions of the segment lengths 

within each composite basis function and are adjusted so that the ordinate to the 
composite triangle is a linear function of displacement along the polygonal line which 
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forms the base of the composite triangle. For example, for five constituent triangles in 
the q th composite triangle of FIG. 20, the weights h q and /z^are 



h q 



(24) 



A* + A* + A* 

1 2 3 



h q 




(25) 



A«+A* + A* 



5 The other weights are computed in a similar fashion. The parameter is the number of 
triangle functions A* employed to represent A . The example composite basis function 

of FIG. 20 is illustrated as the sum of five identical constituent triangles, but, of course, 
the constituents need not be the same if convenience or efficiency dictates otherwise. 
Also, this composite basis function is illustrated without the vector directions associated 
10 with each subdomain. In general the individual straight- wire segments over which a 
composite basis function is defined may each have a different vector direction. Finally, 

the current expanded with a reduced number of unknowns N is 



where A (/) is the q th vector composite basis function defined earlier in (22) and 
15 7 is its unknown current coefficient. It is worth noting that constituent triangles are 



employed above to construct composite triangles but, if desired, they could be used to 
construct other basis functions, e.g., an approximate, composite piecewise sinusoidal 
function. 

If the number of unknowns in a solution procedure is reduced, then, of course, the 
20 number of equations must be reduced too which means that the testing procedure must be 
modified to achieve fewer equations. This is easily accomplished by defining composite 
testing pulses, compatible with the composite basis functions, as a linear combination of 



/(0i(0 = X7 9 A 9 (/)i 9 (0 



(26) 
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appropriately weighted constituent pulses. An example composite test pulse is depicted in 
FIG. 22. Such a p th composite testing pulse is defined by 

N p 

fi^(0 = E<nj(/) (27) 

k=\ 



where the constituent pulses associated with this p th pulse are 

5 n jt (/) = | 1 ' /€(/ /- ?/ ^ (28) 

[0, otherwise 

and shown in FIG. 23. If with every constituent triangle there were associated a 
corresponding constituent pulse, then the testing functions II would overlap, which is 

not desired and can be avoided by selecting the weight w £ to be 0 or 1 depending upon 

whether or not the tf h constituent pulse in II is to be retained. To this end, the inner 

p 

10 product of (13) is modified in the composite testing procedure to become 

(/,n \ = X</:V(0<tf. (29) 

* ' k=\ v 

Now that we have described the new basis and testing functions, we substitute the current 
expansion of (26) into (8) and form the inner product (29) of the resulting expression 

with II for p = 1 , 2,. . ., N . This yields the following matrix equation having a reduced 

p 

15 number (N)of unknowns and equations: 



t z J [/ * ]=[F „ ] (30) 



where 



33 



Ar N q 



4=1 i=l 



— {(A M ' 1 " + A p i " ■ 1 ? )f A 9 0^' 



+ (A " 1 p l 9 + A "\ " \ ") f , A q (l')K(R p )dl' 

^ ic- k- /"+ k+ k+ iV J* / k 



+ ■ 



(31) 



10 



15 



+— tf'KiRpdv 



represents an element of the reduced-rank (n x n) impedance matrix. At this point the 
reader is cautioned to distinguish between the index k which only appears in (31) as a 

subscript and the wave number k = ajjie . The distances R p ± and R p are given in (1 6) 
or (17) with m replaced by index k, and the forcing function is given by 

v =-iuif;w(i[)-\ii)di. 

k=\ v 

One could compute the terms within the reduced-rank impedance matrix directly from 
(31). However, this would require more computation time than needed to fill the original 
impedance matrix of (14) since some constituent triangles within adjacent composite 
basis functions have the same support (FIG. 24). The constituent triangles within the 
overlapping portions of two adjacent composite basis functions differ only in the weight 

h q . Therefore (31) incorporates redundancies which should be avoided. Also, a study of 

(15) and (31) reveals that the term within the braces of (31) is identical to Z mn of (15) if 
subscript i is replaced by n, subscript k by m, and the superscripts p and q are suppressed. 

Hence, the elements Z of the reduced-rank matrix can be computed from the elements 

pq 

Z mn of the original matrix by means of the transformation 



N q 



Z =Yu p Yh q Z pg 

pq f k*-^ i ki 



(32) 
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where Z ^ is a term in the original impedance matrix Z mn of (15). The key to selecting 
appropriate 2 mn term is the combination of indices p, q, k, and /. The index p (q ) 
indicates a group of rows (columns) in J which are ultimately combined by the 

transformation in (32) to form the new matrix. The appropriate matrix element Z^ in 

is determined by intersecting the tf h row within the set of rows identified by index 

p with the i th column of the group of columns specified by index q . Of course the 
groupings of rows and columns are determined when one defines the composite basis and 
testing functions. 



10 A transformation for computing the reduced-rank matrix [Z ] from the 

traditional matrix [Z ] which is more efficient than is the construction of the matrix 

mn 

from (3 1) can be developed. The key to this transformation is (32). First, two auxiliary 
matrices \L 1 and [R 1 are constructed and, then, the desired transformation is 

pm nq 

expressed as 



15 



where 



[Z ] = [L J[ZJ[R J 

pq pm mn nq 



pm 



1 I 1 



2 2 2 
1 2 * 



3 3 3 



p p p 

1 2 N p 



N N N „ 



(33) 



(34) 



and 
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[R ] = 

nq 



h N 2 



1 

h" 

2 



(35) 



It is easy to show that the above matrix transformation is equivalent to (32). 

An alternative development of the transformation, which renders the meaning and 

construction of the matrices [L ] and [R ] more transparent is presented. We begin 

pm nq 

with the traditional NxN system matrix equation, 

[Z J[I n ] = [V], (36) 

which is to be transformed to the NxN reduced-rank matrix equation 



[Z ][I ] = [V ]. 

PQ Q P 



(37) 



The number of unknown current coefficients in the original system of equations (36) is 
reduced by expressing the N coefficients / as linear combinations of the Af coefficients 
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/ (N < N). The / are constructed from the / by means of a scheme which accounts 

n v 7 q n 

for the representation of the composite basis functions in terms of the original triangles 
on the structure. The resulting relationships among the original and the composite 
coefficients are expressed as 

[ / J=[*J^ (38) 

where [7? ] embodies weights of the constituent triangles needed to synthesize 

nq 

composite basis function triangles. The matrix [R ] directly combines unknown current 

coefficients consistent with the composite basis functions to result in a reduced number 
of unknowns. The construction is simple. If the triangle n from the original basis 
10 functions is to be used in the q th composite basis function, the appropriate weight of this 

triangle is placed in row n and column q of [R ] . Otherwise zero is placed in this 

position. Again we point out that a given triangle may appear in more than one composite 
basis function. After substituting (38) into (36) we arrive at a modified system of linear 
equations 

15 [ZJ[*J[7 a ] = FJ (39) 

mn nq q m 

which has a reduced number (N) of unknowns but the original number (N) of equations. 
To reduce the number of equations to N , tested linear equations are selectively added, 
which is accomplished by pre-multiplying (39) by [L ]to arrive at 



20 The identifications, 



and 



[L J[Z J[R J[I ] = [L ][V ]. (40) 

pm mn nq q pnt m 



[Z ] = [L ][Z ][R ] (41) 

pq pm mn nq 



W } = [L ][VJ, (42) 

p pm m 
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in (40) lead to the desired expression (37). The matrix [L ] effectively creates composite 

testing functions from the original testing pulses. If the p th composite testing pulse 
contains the m^ testing pulse from the original formulation, a one is placed in row p and 

column m of [L ] . Otherwise, a zero is placed in this position. 

pm 

5 There are other important considerations in the implementation of this technique. 

Again, we label the number of basis functions in the traditional formulation N and the 
number of composite basis functions N . In the previous section the number of 
constituent triangles for the q'* composite basis function is designated N q Here for ease 

of implementation it is convenient to chose N q to be the same value for every q, which 

10 we designate x ( N q = x for all q). Also, in the present discussion, we restrict x to be one 

of the members of the arithmetic progression 5, 9, 13, 17,...,. With x one of these 
integers, half- width constituent pulses are not required within the composite testing 
functions. TV must be sufficiently large to ensure accurate modeling of the wire geometry 
and vector direction of the current as well as to preserve the numerical accuracy of the 

15 approximations. In addition, N must be large enough to accurately represent the 

variation of the current. A convergence test must be conducted to arrive at acceptable 
values of AT and N . Also, N, N and x must be defined carefully so that a value of x in the 
arithmetic progression will allow an NxN matrix to be reduced to an NxN matrix. 
The following formula is useful for determining relationships between N and N 9 for a 

20 given value of x, in the case of a general three-dimensional curved wire (without 
junctions): 

# = 2— -1. (43) 
r + 1 

For a wire structure with a junction, e.g., a circular loop, where overlapping basis 
functions typically are used in the traditional formulation to satisfy Kirchhoffs current 
25 law, (43) becomes 
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N= -. (44) 

r + 1 

Once N, N and x are determined, it is easy to write a routine which determines the 
original basis and testing functions to be included in the composite functions. This 

information is then stored in the matrices [L 1 and [R ] . 

pm nq 

5 In the above, composite triangle expansion functions are synthesized from generic 

triangle functions but one could as well, if desired, approximate other composite 

expansion functions, e.g., "sine triangles" by adjustment of the coefficients W . Similarly, 

other approximate testing functions could be created by adjustment of the factors u £ . 

Thus, a reduced-rank solution method with composite expansion and testing functions 
10 different from triangles and pulses could be readily created from the techniques discussed 

in this section. Only W and u £ , peculiar to the functions selected in the method to be 

implemented, must be changed in (32) in order to arrive at the appropriate reduced-rank 

matrix elements Z . If [L p J of (34) were replaced by [R J T in (33) where [R J is 

defined in (35) and T denotes transpose, then the resulting reduced-rank matrix [Z ] 

15 would be that for a method which employs composite triangle expansion and 
(approximate) composite triangle testing functions. 
V. RESULTS 

Results obtained by solving the integral equation of (15) with the improved 
solution method developed above are presented in this section as are values of current 

20 determined by the traditional method. In some cases data obtained from the literature are 
displayed for comparison. Results are presented for the wire loop, an Archimedian spiral 
antenna, and several different helical antennas and scatterers. Current values on a small 
wire loop antenna are depicted in FIG. 25. The loop is modeled by 32 linear segments 
(and 32 unknowns) in the traditional solution technique. Also shown are values obtained 

25 from the new solution method with eight composite basis functions (eight unknowns) 
each having five constituent triangles constructed on twenty four linear segments. These 
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current values compare well with those from the traditional solution and with data from 
[2] where the loop is modeled with eight unknowns on quadratic segments. There is 
slight disagreement at the driving point which is to be expected (with eight unknowns) 
near a delta gap source where the current varies markedly. To investigate this discrepancy 

5 we use three triangle basis functions in the vicinity of the delta-gap source and do not 
form composite triangles in this region. The results are shown in FIG. 26. Here the loop 
problem has been solved with 28 unknowns for the traditional method and twelve 
unknowns for the composite basis function solution. It is seen that the agreement is 
excellent even in the vicinity of the delta-gap source. 

10 The improved solution method is applied to a four arm Archimedian spiral 

antenna. This antenna is chosen since it is used in [2] to illustrate the usefulness of the 
quadratic subdomains for wires having significant curvature. A description of the 
geometry of Archimedian spiral antennas is found in [17] and [18]. The antenna is 
excited by a delta gap source on each arm located near the junction of the four arms. The 

15 results presented in this section are for mode 2 excitation [19]. The antenna is also 

modeled by the traditional technique with 725 unknowns on each arm (725*4+3=2903). 
In [17] the authors implement a discrete body of revolution technique so that the number 
of unknowns needed for one arm is sufficient for solving the problem. Since our goal is 
to employ the data of [17] to demonstrate the accuracy of our method and not to create 

20 the best analytical tool for the Archimedian spiral antenna, we solve this problem by 
including the same number of linear segments on each arm and placing overlapping 
triangles at the wire junction to enforce Kirchhoff s current law. In [2] it is found that 
each arm requires 504 linear segments to obtain an accurate solution. They also obtain 
accurate values of the current with 242 quadratic segments. We reproduce these results 

25 with our improved solution method as illustrated in FIGS. 27-29, respectively. The 
number of unknowns for each arm is 725 for the traditional technique and 241 for the 
improved method. In each composite basis function there are five constituent triangles. In 
FIG. 27 the difference in the solution of the current for the two methods is seen to be 
negligible. Good agreement is also achieved for the current magnitude (cf. FIG. 28). A 

30 favorable comparison with data from [2] is observed in FIG. 29. Since the symmetry in 
the geometry is not used to further reduce the number of unknowns required for the 
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structure, the actual number of unknowns in the impedance matrices are 2903 and 967, 
respectively. The computation times for the various routines of the FORTRAN 90 code 
are presented in Table 6 below. All times are for runs on a 375 MHz DEC Alpha 
processor. The time study shows that the reduction technique is successful in 
5 significantly reducing matrix solve time for this four-arm Archimedian spiral antenna. 
A standard linear equation solution method is employed to solve both sets of linear 
equations since the objective of this comparison is to delineate the enhanced efficiency of 
the reduced-rank method. 

TABLE 6 

10 

COMPUTATION TIMES FOR ARCHIMEDIAN SPIRAL 



Event 


Time in Seconds 


Fill matrix N=2903 


1020 


Solve matrix equation N=2903 


1329 


Reduce matrix from 2903 to 967 


5.54 


Solve reduced matrix equation N=967 


45.81 


Traditional method total time 


2349 


Improved method total time 


1071 



Consider next a ten-turn helix having a total wire length of 0.5 A, and illuminated 
by a plane wave. The geometry of the helical scatterer is depicted in FIG. 30. The current 

15 shown in FIG. 3 1 is "converged" when the number of unknowns in the traditional 

solution technique reaches 259. Thus one concludes that 260 linear segments are required 
to accurately represent the geometry of this structure and vector nature of the current. We 
determine convergence by examining the real and imaginary parts of the current along 
the structure. When changes in the current are sufficiently small as the number of 

20 segments is increased, convergence is assumed [2]. The results of a convergence test 
show that an accurate solution of the current can be achieved with 51 composite basis 
functions. The number of constituent triangles in each basis function in this case is nine. 
We note that the solution with 27 composite basis functions differs only slightly from the 
converged solution. 

25 The current is shown in FIG. 32 for another helical scatterer of geometry similar 

to that described above and subject to the same excitation and geometry similar to that 
described above. The circumference of each turn of this ten-turn helix is 0.035X making 
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the total wire length 035X. These results are given as an example to illustrate that the 
composite basis function scheme works well with curved- wire structures having a wire 
length which is not an integer multiple of half wavelength. 

The data of FIG. 33 are for a 50-turn helix having a total wire length of 2X and 

5 illuminated by a plane wave traveling in the positive x direction. One sees that 27 
unknowns are adequate to accurately represent the current along the helix. However, 
1483 unknowns are required in the traditional solution method since many linear 
segments are required to define the 50-turn structure and the vector properties of the 
current. In this example there are 105 constituent triangles in each composite basis 

10 function. Table 7 below shows the computational savings enjoyed by the method of this 
invention. 

TABLE 7 

COMPUTATION TIMES FOR FIFTY-TURN 
15 HELIX 



Event 


Time in Seconds 


Fill matrix N=l 483 


300 


Solve matrix equation N=1483 


267 


Reduce matrix rank from 1483 to 27 


1.84 


Solve reduced matrix equation N=27 


Negligible 


Traditional method total time 


567 


Improved method total time 


302 



Next we illustrate the prowess of the solution technique for helical antennas. 
Specifically the data presented in FIG. 34 and FIG. 35 are for helical antennas driven 
above a ground plane by a delta gap source. The geometry of the helix is given in FIG. 30 

20 and the ground plane is located at z = 0. The data of the improved method compare well 
with those of the traditional solution technique, but, again, there is a slight difference in 
the currents at the ground plane due to the nature of the delta gap source. In each of these 
figures the number of unknowns given is the number for the structure plus its image, but 
' data are plotted only for the part of the structure above the ground plane. Since there are 

25 many turns, the number of segments needed to represent the geometry of the antenna and 
its image is large. The number of unknowns is reduced from N=917 in the traditional 
method to N=53 in the improved technique. Of course, one could employ image theory to 
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modify the integral equation which could be solved by the new method with an even 
more dramatic savings in computer resources. 

The last example is a five-turn helical antenna over an infinite ground plane, 
driven by a delta gap source. This structure is included here because it is used in [6] to 
5 exhibit the accuracy of a technique employing basis and testing functions defined over 
arcs of constant radii. It is modeled by straight wire segments in [20]. In [6] the authors 
discretize the antenna into fifteen arcs and then compare solutions of 135 unknowns with 
forty-five unknowns. They find that forty-five unknowns is enough to obtain an accurate 
solution for the current when the geometry is defined by arcs. We reproduce these results 

10 except that the antenna geometry is defined by many straight wire segments. In the 

method of this invention we include the unknowns on the image (269 unknowns on the 
antenna plus its image corresponds to 135 unknowns on the antenna above the ground 
plane). Likewise, 89 unknowns on the antenna and image are equivalent to 45 unknowns 
on the antenna. We find that helical antennas require a minimum of 25 unknowns per turn 

15 in the traditional solution technique in order to represent the geometry. In order to reduce 
the number of unknowns over the antenna and its image from 269 to 89, each composite 
basis function is constructed with 5 constituent triangles. A qualitative comparison of our 
data and that of [6] suggests agreement in the two methods. 

20 VI. CONCLUSIONS 

The solution method presented in this specification is very simple and practicable 
for reducing the rank of the impedance matrix for curved-wire structures. It should be 
mentioned that rank reduction is realized only when the number of segments needed to 
model the geometry and vector direction of the current exceeds the number of unknown 

25 current coefficients necessary to characterize the variation of the current. We define 

composite basis and testing functions as the sum of constituents over linear segments on a 
wire and arrive at a new impedance matrix of reduced rank. It is shown how this reduced- 
rank matrix can be determined from the original impedance matrix by a matrix 
transformation. Thus one advantage of this technique is that it can be applied to almost 

30 any existing curved- wire codes which define basis and testing functions over straight- 
wire segments or curved- wire segments. 
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Dramatic savings in matrix solve time are realized for the cases of the four-arm 
Archimedian spiral antenna and the helical antenna. The benefits for reducing unknowns 
on, for example, a helical antenna become much more significant as the number of turns 
increases. It should be pointed out that this method does not reduce matrix fill time since 

5 the elements of the original impedance matrix are computed as a step in the determination 
the elements of the reduced-rank matrix. Problems involving large curved-wire structures 
can be solved readily by this method, e.g., a straight wire antenna loaded with multiple, 
tightly wound helical coils and an array of Archimedian spiral antennas. The principles 
described here can be used in addition to other methods such as those based upon 

10 iteration. 
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An exemplary genetic algorithm that can be used in accordance with the subject 
invention to obtain optimal antenna parameters for given design criteria is included in the 
computer program listing appendix provided on compact disc and is incorporated by 
reference herein. 



The following portion of the specification, especially with reference to FIGS. 37- 
47C, respectively, particularly concerns bandwidth enhanced normal mode helical 
antennas. It begins by setting forth the objectives, considerations, and questions 
addressed in the beginning stages of development of the present invention. The effects of 

20 different physical antenna parameters on antenna performance are addressed by showing 
the effect in the VSWR by these variations. 

The remainder of the discussion with respect to FIGS. 37-47C, respectively, 
shows several different antenna designs and in graphical form illustrates the respective 
performance of each. A straight wire antenna, a simple helix, and a triple helix are all 

25 examined. Each antenna is modified by the addition of various parasitic elements. The 
characteristics of each of these antennas are then illustrated. The VSWR, directivity, and 
input impedance are shown so that the different antennas having different combinations 
of parasitic elements can be analyzed effectively. Results obtained from the different 
antenna combinations are then summarized and conclusions drawn from these results are 

30 set forth. Such results illustrate the initial indications that bandwidth improvements 
could be made by the addition of these parasitic elements. 

Objectives for the subject antenna include it being low-profile, omnidirectional and 
broadband. Design considerations include: (a) the helix can be made shorter by adjusting 
the pitch, (b) normal mode helix has narrow bandwidth, and (c) parasitic elements 
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increase the bandwidth of straight wires. Questions addressed in the subject analysis 
include determination of whether the bandwidth of the normal mode helix can be 
improved with parasitic elements and if so, what are suitable structures for the parasites. 
FIG. 37 presents numerical results of the effect of pitch angle on helix VSWR for 

5 a helical antenna having a total wire length of 75cm and a wire radius of 0.5cm. The 
different curves plotted in the graph of FIG. 37 are for antennas having a helix pitch (0) 
of 30, 40, 50, 60 and 70 degrees. For an antenna having a pitch of 30 degrees, a height 
reduction of 45.0% and a total height of 41.3cm is achievable. For an antenna having a 
pitch of 40 degrees, a height reduction of 32.0% and a total height of 50.9cm is 

10 achievable. For an antenna having a pitch of 50 degrees, a height reduction of 21 .0% and 
a total height of 59.2cm is achievable. For an antenna having a pitch of 60 degrees, a 
height reduction of 12.0% and a total height of 66.0cm is achievable. For an antenna 
having a pitch of 70 degrees, a height reduction of 5.4% and a total height of 70.9cm is 
achievable. 

15 FIG. 38 presents numerical results of the effect of wire radius on helix VSWR for 

a helical antenna having varied diameters, including 0.4 cm, 1 cm and 2 cm, as depicted 
in the graph. The helix geometry of such antennas are characterized by a height of the 
straight wire base being 7.5 cm, the circumference of one turn in the helix being 15cm, 
total number of turns being 4.5 with a pitch angle of 40 degrees, and a total wire length of 

20 75cm. 

FIG. 39 provides a block diagram representing exemplary steps in a procedure for 
efficient optimization of a helix with parasitic elements. The evaluation steps are done 
for each antenna in the sample population. In the steps of FIG. 39, Z mn is the reduced- 
rank impedance matrix used for curved wires. 

25 FIG. 40A provides numerical results comparing the VSWR versus frequency for 

an open sleeve monopole antenna, such as depicted in FIG. 40B, and a regular straight 
wire antenna. For a single wire antenna, a VSWR less than 3.5 is achieved on a 
frequency range from 85 MHz to 1 12 MHz for a bandwidth ratio of 1 .32: 1 . For an open 
sleeve monopole antenna, a VSWR less than 3.5 is achieved on a frequency range from 

30 90 MHz to 1 72 MHz for a bandwidth ratio of 1 .9: 1 . 
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FIG. 41 A provides numerical results comparing the VSWR versus frequency for a 
straight wire antenna having four parasites, such as depicted in FIG. 4 IB, and a regular 
straight wire antenna. For a single wire antenna, a VSWR less than 3.5 is achieved on a 
frequency range from 85 MHz to 1 12 MHz for a bandwidth ratio of 1 .32: 1 . For a straight 
5 wire antenna having four parasites, a VSWR less than 3.5 is achieved on a frequency 
range from 90 MHz to 1 85 MHz for a bandwidth ratio of 2.05: 1 . 

FIGS. 42A through 45B present numerical results for various helix antenna 
embodiments. The basic geometry of the helical antenna is the same as previously 
described for the antennas modeled in FIG. 38 and having a wire diameter of 1 cm. 
10 FIG. 42 A provides numerical results comparing the VSWR versus frequency for a 

helix antenna with two straight wire parasites, such as depicted in FIG. 42B, and a regular 
helix antenna. For the helix antenna having two parasites, a VSWR less than 3.5 is 
achieved on a frequency range from 1 12 MHz to 208 MHz for a bandwidth ratio of 
1.86:1. 

15 FIG. 43 A provides numerical results comparing the VSWR versus frequency for a 

helix antenna with four straight wire parasites, such as depicted in FIG. 43B, and a 
regular helix antenna. For the helix antenna having four parasites, a VSWR less than 3.5 
is achieved on a frequency range from 1 12 MHz to 250 MHz for a bandwidth ratio of 
2.23:1. FIG. 43C provides a graph of the directivity versus frequency for antennas with 

20 varied pitch angles (0) for <|)=0. FIG. 43D provide a graph of the directivity in the H- 
plane versus § when f = 190 MHz and 9 = 90 degrees. 

FIG. 44A provides numerical results comparing the VSWR versus frequency for a 
helix antenna with two helical parasites, such as depicted in FIG. 44B, and a regular helix 
antenna. For the helix antenna having two helical parasites, a VSWR less than 3.5 is 

25 achieved on a frequency range from 1 12 MHz to 208 MHz for a bandwidth ratio of 
1 .86: 1 . Real and imaginary components for the input impedance of the antenna 
represented in FIG. 44B is displayed in FIG. 44C, and FIG. 44D charts the directivity 
versus frequency for antennas having different pitch angles (6). 

FIG. 45 A provides numerical results comparing the VSWR versus frequency for a 

30 helix antenna having respective inner and outer helical parasites, such as depicted in FIG. 
45B, and a regular helix antenna. For the sleeve helix antenna having inner and outer 
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parasites, a VSWR less than 3.5 is achieved on a frequency range from 101 MHz to 182.5 

MHz for a bandwidth ratio of 1 .8 1 : 1 . 

FIGS. 46C through 47D, respectively, present numerical results for various triple 

helix antenna embodiments, such as represented in FIGS. 46A and 46B. The basic 
5 geometry of the triple helix antenna is the same as previously described for the antennas 

modeled in FIG. 38 and having a wire diameter of 1 cm, except the triple helix has 13.5 

turns (4.5 for each helix). 

FIG. 46C provides numerical results comparing the VSWR versus frequency for a 

triple helix antenna and a single helix antenna. FIG. 47A provides numerical results for a 
10 triple helix antenna having four straight wire parasites, such as depicted in FIG. 47B, a 

triple helix antenna and a single helix antenna. For the triple helix antenna having four 

parasites, a VSWR less than 3.5 is achieved on a frequency range from 110 MHz to 380 

MHz for a bandwidth ratio of 3 .45 : 1 . FIG. 47C provides a graphical representation of 

antenna directivity versus frequency for the triple helix antenna with different values for 
15 the helix pitch angles (0). FIG. 47D provides a graphical representation of antenna 

directivity versus frequency for the triple helix antenna with parasites (such as in FIG. 

47B) with different values for the helix pitch angles (0). 

A summary of the results determined from the numerical data provided in FIGS. 

40A through 47D, respectively, is now presented in Table 8 below. These results show 
20 that parasitic straight wires and helices are useful for improving the bandwidth of helical 

antennas. Also, the triple helix has reduced VSWR over the frequency band which 

makes the structure more amenable to improvement by parasites. 



Driven Element 


Number of 
Parasites 


Types of Parasites 


Bandwidth Ratio 


Straight Wire 


2 


Straight wire 


1.90:1 


Straight Wire 


4 


Straight wire 


2.05:1 


Helix 


2 


Straight wire 


1.86:1 


Helix 


4 


Straight wire 


2.23:1 


Helix 


2 


Helix (same cylinder) 


1.86:1 


Helix 


2 


Helix (different cylinders) 


1.81:1 


Triple Helix 


4 


Straight wire 


3.45:1 



Table 8 

25 
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ABSTRACT OF THE DISCLOSURE 



A method for applying an algorithm to facilitate the design of wideband 
omnidirectional antennas, and the design of sleeve cage monopole and sleeve helix units 

5 includes rapid resolution of a complex relationship among antenna components to yield 
an optimal system. A genetic algorithm is used with fitness values for design factors 
expressed in terms to yield optimum combinations. Cage antennas are optimized via a 
genetic algorithm for operation over a wide band with low VSWR. Genetic algorithms 
and an integral equation solver are employed to determine the position and lengths of 

10 parasitic wires around a cage antenna in order to minimize VSWR over a band. The cage 
may be replaced by a normal mode quadrifilar helix for height reduction and with re- 
optimized parasites. 



